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ABSTRACT 

The mass function of galaxy clusters is a powerful tool to constrain cosmological param- 
eters, e.g., the mass fluctuation on the scale of 8 h^ 1 Mpc, erg, and the abundance of total 
matter, Vt m . We first determine the scaling relations between cluster mass and cluster rich- 
ness, summed r-band luminosity and the global galaxy number within a cluster radius. These 
relations are then used to two complete volume-limited rich cluster samples which we ob- 
tained from the Sloan Digital Sky Survey (SDSS). We estimate the masses of these clusters 
and determine the cluster mass function. Fitting the data with a theoretical expression, we 
get the cosmological parameter constraints in the form of a%{Vl m / '0 .3) Q = (3 and find out 
the parameters of a =0.40-0.50 and f3 =0.8-0.9, so that a 8 =0.8-0.9 if Q. m = 0.3. Our cr 8 
value is slightly higher than recent estimates from the mass function of X-ray clusters and the 
Wilkinson Microwave Anisotropy Probe (WMAP) data, but consistent with the weak lensing 
statistics. 
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1 INTRODUCTION 

Precise determination of cosmological parameters is an impor- 
tant goal in astrophysics. In the linear theory, the present root- 
mean-square (rms) mass fluctuation on the s cale of 8 ft" 1 Mpc , 
erg, is one of fundamental parameters (see ISpergel et alj 120031) 
to describe the power spectrum of mass fluctuations in the uni- 
verse. It is one of key parameters in the large scale struc- 
ture simulations (e.g., Ijenkins et alj|l998l) . The ag can be deter- 
mined by galax y-galaxy correlations (e.g.. iTegmark et all 120041 : 
ICole et al.l2 005). fluct uation s in the cosmic microw ave background 



(Sper gel et alj l200j 120071: iKomatsu et alj |2009l). gravitational 



lensing statistics (e.g ., iHoekstra et alj|2006l ; iKit ching et al.ll2007l; 



Benjamin et alj l2007h. cluster mass function (e.g., White et all 



19931: iBahcall & Farll998l;lReiprich & Bohringej|2002h . Lyq for 
est jjena et al. 2005 : McDon ald et alj20050 and galaxy peculiar ve- 
locities ( Feldm an et alfeoolf) ! 

The cluster mass function can be determined by the estimated 
masses for a sample of clusters (e.g., Dahle 2006), or by the X-ray 
luminosity and temperature function wit h a prior scaling relation 
dViana & Lidd5ll996l : lAllen et alj|2003t) . Fitting the cluster mass 
function with a theoretical expression can provide constraint on a$. 
Generally, erg is coupled with Q. m , the abundance of present total 
matter, in the form of crg(f2 m /0.3) O! = /3. Previous studies have 
found a in the range 0.3-0.6 and /? in the range 0.6-1 .2 (see Table|2] 
in Section[4]l. The determined erg in recent years (2002-2009) has a 
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mean value of 0.73±0.05 as suming Q, m = 0.3, w hich is in agree- 
ment with the WMAP data (Komatsu et al. 200 9j), but lo wer than 
those by weak lensing statistics ( Hetterscheidt et al. 2007), galaxy- 
galaxy correlations (Tegmark et al. 2004; Cole et alj2005h and Lya 
forest Jjena et alj|2005l ; lMcDonald et alj2005h . 



The amplitude of cluster mass function has large uncertain- 
ties, mainly cause d by the uncertain normalization of the mass scal- 
ing relation (e.g., lHenrv||2004h . Other uncertainties come from the 
scatter of mass scaling relation and the incompleteness of t he X- 
ray flux-limited cluster samples (Reiprich & Bohringer 2002]). The 
cluster m ass function may be und erestimated if only X-ray clusters 
are used. lErben et alj (2000) and iDahle et alj d2003l) have noticed 
the existence of a cla ss of X-ray-underluminous massive clusters. 
IPopesso et alj d2007al) found that 40% of Abell clusters have a low 
level or no detection in X-rays. A large complete volume-limited 
sample of clusters is crucial for the purpose. Using the photomet- 
ric redshifts of galaxie s, we found 39, 668 clusters in the redshift 
range 0.05 < z < 0.6 dWen et alj|2009h . Clusters are approximate 
volume-limited complete in the redshift range 0.05 < z < 0.42. 
The richnesses and the summed luminosities of clusters are esti- 
mated from their luminous members, and they are tightly related to 
cluster mass. In Section 2, we carefully determine the scaling rela- 
tion for cluster mass. In Section 3, we get the cluster mass function 
for a local sample of clusters and a sample at mediate redshifts, and 
then fit the cluster mass function with a theoretical expression for 
constraints on cosmological parameters, Q. m and <rg. Discussions 
and conclusions are given in Section 4. 
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Table 1. Cluster masses from literature and the mass tracer values for 53 clusters (richness R > 8 and 0.03 < z < 0.3, sorted with z) in the 
field of the SDSSDR6. 



Name 


RA. 


Decl. 


2 


R 




GGN/r GGN 


Mvir 




(deg) 


(deg) 






(1U ft Lq) 


(Mpc ) 


/ 1 nl4 i, — 1 7i/f \ 
(_ 1U il M0) 


Abell 2199 


247.15930 


39.55121 


0.030 


17.76 


28.49 


18.10 


o oq+0.38 

z - oa -0.38 

3 42+ ' 26 
°-^ z _0 26 
















Abell 2052 


229.18536 


7.02162 


0.035 


9.50 


21.97 


11.62 


I 4V+0.28 
-0.28 
i co+0.05 
i - &S -0.06 
















Abell 2063 


230.77209 


8.60922 


0.035 


15.30 


22.96 


19.93 


9 01+O.I8 
z,ol -0.16 

9 70+O.42 


Abell 2147 


240.57094 


15.97465 


0.035 


12.22 


22.50 


14.67 


2 46+ ' 83 
z -™-0.52 


Abell 215 lw 


241.14914 


17.72156 


0.037 


11.08 


22.16 


13.29 


-, 90+O.O9 
1 ' zo -0.09 


MKW9 


233.13339 


4.68100 


0.040 


8.88 


15.77 


11.44 


1 01+ ' 25 

1 06+ ' 27 
±.uo_ 27 
















Abell 1983 


223.23048 


16.70286 


0.044 


9.27 


11.76 


11.45 


l-39toT 3 
l-41±g : gg 


Abell 160 


18.24822 


15.49129 


0.045 


16.61 


34.96 


16.54 


0.96l° : « 


Abell 85 


10.46029 


-9.30312 


0.056 


23.61 


54.65 


21.85 


4.23+^ 

5 55+O.58 
—0 53 
















Abell 1991 


223.63122 


18.64232 


0.059 


19.22 


39.13 


17.48 


1 05+0:19 

1 44+O.17 
-"-•^-O.IT 
1 59+O.21 
1 ' oz -0.21 

1 66+ ' 32 
±.uo_ 24 
















Abell 1795 


207.21877 


26.59293 


0.063 


19.56 


39.70 


19.56 


5 49+0.65 
°-^ z -0.65 

7.58 1 7q 
7.86+» : 7 ° 


Abell 2092 


233.31403 


31.14515 


0.067 


9.27 


15.93 


11.45 


1 09+ ' 23 


Abell 2065 


230.60008 


27.71436 


0.072 


40.86 


76.31 


40.86 


12 ii+ 15 -48 
iz.ii_ 4 gg 


ZwCl 1215 


184.42134 


3.65584 


0.075 


20.69 


43.81 


18.71 


7.54+fH 


Abell 1800 


207.34822 


28.10732 


0.075 


19.77 


44.73 


18.83 


4.14+^ 


Abell 1775 


205.45477 


26.37347 


0.076 


14.52 


32.62 


17.76 


n7 +0.43 

°- ul -0.29 


Abell 2029 


227.73376 


5.74478 


0.077 


24.87 


60.63 


27.47 


7 06+ ' 69 
' - UO _0.54 
30+1-93 

10-55±i;81 
















Abell 2255 


258.11996 


64.06072 


0.080 


40.23 


83.40 


33.51 


7 +1.05 


Abell 1650 


194.67288 


-1.76146 


0.084 


17.08 


42.92 


20.68 


8 1 q+4-21 
°- id -2.53 


Abell 1692 


198.05661 


-0.97448 


0.085 


14.40 


26.42 


14.73 


I ig+0.38 


Abell 1750 


202.71080 


-1.86197 


0.088 


24.35 


58.05 


27.58 


4 78+ 2 - 64 
lo -2.64 


Abell 2142 


239.58334 


27.23341 


0.090 


39.33 


80.90 


39.88 


ll.OOl 2 '^ 

11.95+1 24 


Abell 2244 


255.67705 


34.05999 


0.097 


31.22 


62.35 


28.67 




Abell 2034 


227.54883 


33.48646 


0.113 


32.00 


68.97 


31.08 


7 -| 7 +4.3{l 

— 4 30 


Abell 1068 


160.18541 


39.95313 


0.138 


16.38 


38.87 


16.38 


4 87+° :42 

^•°'-0.B3 

5 19+O.45 
°- ±z -0.45 
















Abell 1413 


178.82501 


23.40491 


0.143 


37.17 


93.39 


41.82 


5 70+0.57 
°'' U -0.53 
5 07+O.59 
O-<5'_0.59 

8 -' 2 -0.70 
76 +2.54 
a -'°-1.89 
05+1.04 
















RXJ1720. 1+2637 


260.04184 


26.62557 


0.164 


26.87 


69.11 


27.30 


c 00+4. 98 
5.88_ 4 60 

5 84+ 4 ' 30 
°'°^-4.30 
















Abell 1914 


216.48611 


37.81645 


0.171 


34.16 


84.15 


38.07 


4 fifi +3 ' 78 
4,DD -3.49 

18.40l 2 ; 2 ° 
















MS 0906.5+1110 


137.30312 


10.97475 


0.176 


36.00 


119.08 


38.61 


8.30+j; 30 , 



Method Ref. 



X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

WL 

X-ray 

WL 

X-ray 

WL 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

X-ray 

WL 

WL 

WL 

WL 

WL 

X-ray 

WL 



3 
4 
3 
4 
2 
7 
6 
9 

10 
9 
8 

10 
2 

11 
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Table 1. - continued 



Name 



R.A. 
(deg) 



Decl. 
(deg) 



R L r GGN/r GGN 

(10 10 h- 2 L Q ) (Mpc- 1 ) 



( 10 



M vir 
14 h~ 



Method Ref. 



1 M Q ) 



Abell 1689 



Abell 963 



Abell 1682 
Abell 1763 

Abell 2219 

Abell 2111 
Abell 267 



Abell 1835 



197.87291 -1.34108 0.184 62.74 



154.26515 39.04705 0.206 43.00 



Abell 1423 179.32219 
RXJ1504. 1-0248 226.03130 
Abell 773 139.47261 



196.70833 
203.83372 



234.91872 
28.17483 



Zw 1231.4+1007 188.57277 
MS 1231.3+1542 188.48055 
RXJ2129.6+0005 322.41649 



210.25863 



33.61092 0.214 20.78 
-2.80460 0.215 20.08 
51.72704 0.217 66.83 



46.55927 0.226 40.76 
41.00115 0.228 34.10 



250.08253 46.71148 0.228 45.62 



34.42426 0.229 47.58 
1.00711 0.230 31.53 



9.76623 0.231 32.00 
15.43305 0.234 13.82 
0.08920 0.234 24.00 



2.87846 0.252 51.57 



MS 1455.0+2232 224.31295 22.34288 0.258 26.00 



Abell 2631 354.41554 0.27138 0.277 46.22 

Abell 1758N 203.16007 50.55992 0.279 42.24 



136.85 



93.83 



64.82 
73.32 
138.50 



112.14 
109.20 

94.82 

112.27 
86.85 



98.28 
35.53 
80.65 



129.59 



72.03 



113.32 
107.00 



50.32 


10.50tJ|J 


WL 


12 




1 1 fiQ +1 ' 80 


SL+WL 


13 




io 24+3.99 
10.z^_ 3 gg 


X-ray 


14 




13.431°$ 


X-ray 


2 




13.5ltl;« 


WL 


15 




14 70+ 1 - 40 


WL 


16 




17 KK+3.04 
l/.00_ 3 04 


WL 


17 




20.53ti^ 


WL 


18 




25.50+i;3_°, 


WL 


11 


38.52 




WL 


17 




fi Kn+ 2 ' 00 
D.OU_ 1 gg 


WL 


11 






X-ray 


14 




7 on+ 2 - 45 

' - uu -1.60 


X-ray 


7 




8-01_ 6 ' 36 


WL 


9 




a oq+7.46 
8-3»_7.46 


WL 


10 


21.16 


14.98l 7 ;f 7 


WL 


10 


26.82 


15.101^° 


X-ray 


7 


55.08 


10.63±1;1? 


X-ray 


14 




1 ^ QQ+6-95 
1J -°°-6 .95 


WL 


10 




94 7f)+ 9 - 44 


WL 


9 


39.49 




WL 


9 




5.49l| : ^ 


WL 


10 


32.77 


6.23±i;l| 


X-ray 


14 




s.9itif 


WL 


9 




10.54±g;i« 


WL 


10 




12.28 + 2'37 


WL 


17 




13.50ti-;™ 


WL 


11 


45.62 


8.19lg:S 


WL 


10 




°-°'-5.4S 


WL 


9 




11.30±1?° 


WL 


11 




18.66l^ 4 


WL 


17 


43.38 


6.12±1;|* 


WL 


10 




6-92ti;i 


WL 


9 


32.62 


2 3fi +1 ' 29 
.29 


WL 


17 






X-ray 


14 




(.OU_ 2 2 o 


WL 


11 




14.98±|-|| 


WL 


10 




16.34l|; 7 7 


WL 


9 


36.99 


Q 1 ir + 2.65 
°- 1J -2.65 


WL 


10 


17.40 


i Kn+0.90 
± - JU -0.90 


WL 


11 


20.45 


o.yo_ 2 70 


X-ray 


7 




94 

7-5i±l;lt 


WL 


10 




15.491^1 


WL 


9 


47.96 


8.19l 3 3 ;° 7 


WL 


12 




8.41+ 2 ' 5 57 


X-ray 


14 




10.61±|;«® 


WL 


10 




15.62ll;l« 


WL 


9 




17.00tii° 


X-ray 


7 




24.21±ij;™ 


WL 


17 


26.00 




WL 


9 




7 30+1-90 
' - ou -1.80 


WL 


11 




12.88t«$ 
6.12±1|| 


WL 


10 


44.94 


WL 


10 


41.36 


c 2 fi+ 5 - 70 
°- zo -5.70 


WL 


8 




26.9lt|;g 


WL 


10 




39.09ti 2 ;^ 


WL 


9 
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Table 1. - continued 



Name 


RA. 

(deg) 


Decl. 
(deg) 


2 


R 


(10 10 h- 2 Lq) 


GGN/r GGN 
(Mpc- 1 ) 


M vir 
( 10 14 h- 1 Mq) 


Method 


Ref. 


Abell 697 


130.73982 


36.36646 


0.282 


27.72 


77.66 


37.83 


22 Q9+ 9 - 78 

26 10+ 9 - 54 
zo.±u_ 9 54 


WL 
WL 


9 
10 


Abell 959 


154.39984 


59.56710 


0.285 


52.77 


139.82 


48.82 


17.38l 6 5 ;« 


WL 


9 


Abell 611 


120.23674 


36.05655 


0.288 


27.86 


91.23 


30.92 


6.18±?;|» 
6 54+ 4 ' 43 


X-ray 


7 
















WL 


10 
















6.90tHl 


WL 


9 


Abell 781 


140.20117 


30.47176 


0.288 


27.17 


77.19 


33.22 


LZ -°< -5.86 


WL 


10 


Zw3146 


155.91515 


4.18629 


0.291 


31.76 


70.81 


30.33 


n si+ 5 - 44 

ll.dl_ 5 44 
13.99+™? 


WL 
WL 


10 

9 


Zw 1459.4+4240 


225.34604 


42.34448 


0.292 


27.67 


81.26 


29.89 


in KO+6-00 

iu.oy_ 6 00 


WL 


10 


Abell 1576 


189.24684 


63.18658 


0.300 


42.00 


156.45 


38.40 


I6.0l±«;|g 


WL 
WL 


9 
10 



Note for Method: X-ray stands for the mass determined by X-ray measurements; W L st ands for weak lensing, WL+SL sta nds for weak lensing 
combined with strong lensing; References for mass estimates: [1] Xu et al. 12001), [2] Reiprich & Bohringer (2002), [3] Pointecouteau et al. 
120051). 14llArnaud et alj l2005t) 15llStin et alj <2009l). r6llvikhlinin et alj fcOodTmlSchmidt & AllerJ 120071). 18llOkabe & Umetsul 120081). 191 
iPedersen & Dahlel i200T), 1101 iDahld J2006l), iTll iHoekstr j fe007l). ri 2l TCorless et alj <2009l). 11 31 iLimousin et alj <2007l) 1141 IZhang eTall 
J2007I) . 115l lBroadhurst etal] j2005t) . ri6l lUmetsu & Broadhurstl )2008l) . 117l lBardeau et alj J2007l) . 118l lHalkola et aljj2006l) . 




R L r (10 10 h" 2 L ) GGN/r GGN (Mpc" 1 ) 

Figure 1. Correlations between cluster mass Af v j r and richness R, summed luminosity L r and GGN/tqqpj for 24 nearby clusters (R ^ 8 and z < 0.1). 
The solid line is the best fit as given in Equation (D-^3j. 



Throughout this paper, we assume a flat ACDM cosmology, 
taking H =100/ikms _1 Mpc" 1 , with h = 0.72, Q m = 



2 MASS SCALING RELATIONS FOR CLUSTERS 

We identified 39,668 clusters from the SDSS DR6 by discrimina- 
tion o f luminous member galaxies with following steps i Wen et all 
2009). First, we assume that each galaxy at a given photometric red- 
shift z is the central galaxy of a cluster candidate, and we count the 
number of luminous "member galaxies" of M r ?C —21 within a ra- 
dius of 0.5 Mpc and a photometric redshift gap of z ± 0.04(1 + z). 
We set Az = 0.04(1 + z) for the gap to allow variable uncer- 
tainties of photometric redshifts at different redshifts. Second, we 



define the center of a cluster candidate to be the position of the 
galaxy with a maximum number count. The cluster redshift is es- 
timated to be the median value of the photometric redshifts of the 
recognized "members". Third, for each cluster candidate at z, all 
galaxies within 1 Mpc from the cluster center and z ± 0.04(1 + z) 
are assumed to be the member galaxies. Their absolute magnitudes 
are re-calculated with the cluster redshift. Finally, a cluster at z is 
identified when the number of member galaxies of M r ^ —21 
reaches 8 within a projected radius of 0.5 Mpc and z ± Az. Monte- 
Carlo simulations show that the detection rate is more than 90% for 
massive clusters (richness R ^ 16.7) if the redshift uncertainty of 
cluster galaxies is about 0.03(1 + 2). 

We defined the cluster richness, R, to be the total number of 
galaxies(Af r < —21) within a radius of 1 Mpc and z±0.04(l + z) 
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Figure 2. Correlations between cluster mass Af v i r and richness R, summed luminosity L r and GGN /tggn for 17 clusters in the redshift range 0.17 < 
z < 0.26. The black dots are the clusters with more than three estimates of their masses. The solid lines are the same shown in FigurefT] The dashed line is 
the new scaling relation with the same slope but different offsets determined from the data. 



after subtracting the local background, i.e., the average number of 
luminous galaxies. The summed r-band luminosity of each cluster, 
L r , is calculated as the total luminosity of member galaxies within 
the region also after subtracting the background. From the radial 
distribution of member galaxies, we got the cluster radius, tggn * 
where the density of galaxies is as low as background. Here, we 
defined the Gross Galaxy Number (GGN) of a cluster as the to- 
tal number of luminous galaxies (M r ^ —21) within the radius 
tggn and the redshift gap of z ± 0.04(1 + z) after subtracting 
the local background. It has been known for a long time that the 
cluster richness and summed luminosity ar e related to cluster mass 
dGirardi et all 120021 ; IPopesso et al.ll2007bl) . hence they can be the 
tracers of cluster mass. The GGN/rGGN is related to the ampli- 
tude of cluster-galaxy cross- correlation since the cor relation is de- 
scribed by £(r) oc r -2 (e.g jLilie & Efstathioul 19881) . We find that 
GGN/rGGN can also be the tracer of cluster richness. 

Cluster mass can be determ ined by the velocity dispersion of 
member galaxies dZwickvl 19331) . However, velocity measurements 
can be corrupted by projection effects that might be difficult to di- 
minish in practice. The err or on the individual measur ements can 
introduce a significant bias dvon der Linden et alj2007h . Under the 
assumption of hydrostatic equilibrium, the X-ray method can deter- 
mine mass distribution of a cluster to a large radius. The assump- 
tion is invalid for cl usters with subs tructures, inducing an under- 
estimation of mass dSchindleilll996h . Weak gravitational lensing 
recently becomes a sophisticated method to estimate cluster mass 
without assumptions on dynamical state of a cluster. The uncer- 
tainty of mass mainly comes from the difficulty in measuring the 
image distortions of the faint background sources. We collect the 
cluster masses estimated by X-ray and weak lensing methods from 
literature (see Table[T). Usually, cluster masses are denoted as Ma 
which is the mass within a radius r& interior to which the mean 
density is A times the critical density of the universe. For cosmol- 
ogy with Q m — 0.3, the virial mass is calcula ted within the ra- 
dius t a, here A = 101, so that M v i r = Mioi dKitavama & Sutol 
1996). Pre vious studies usually provided the mass within r2oo o r 
r 50 o (e.g.. iReiprich & B6hringerll2002l ; IPedersen & DahleH2007t) . 



Here, we convert t he mass of M2 and M500 to the virial mass 
Mvir according to IShimizu et alj d2003l) . We will discuss later 
the influence on our result from a possible bias conversion. For 
each cluster with mass estimated, we calculate the cluster rich- 
ness, the sum med r-band lumin osity and GGN/rGGN following 
the method of IWen et alj d2009h . Only clusters of richness R ^ 8 
are listed in Table [JJ since the uncertainties of R and the summed 
luminosities become larger for clusters with a smaller R. 

We notice that clusters with estimated masses preferentially 
have low (2 < 0.1) and mediate (~ 0.2 < 2 < 0.25) redshifts (see 
Table[JJl. To minimize the uncertainty, we determine the scaling re- 
lations between the masses and observational tracers for clusters 
in the two small redshift ranges independently. This is because the 
discrimination of member galaxies (e.g., completeness or contami- 
nation rate) may be different for clusters at different redshifts, and 
the systematic bias can be ignored in such a small range. In the low 
redshift range (z < 0.1), the masses of many clusters are available 
and distributed in a large mass range, which is good for determina- 
tion of the scaling relations. We get 15 clusters of 0.05 < z < 0.1. 
We also include 8 clusters of 0.03 < 2 < 0.05 and one cluster of 
2 = 0.113 to derive the scaling relations at the low redshift range. 
Several clusters have multiple estimates for mass from literature, 
we adopt the median value or the average of two middle ones for 
even measurements. 

The mass-richness relation, i.e., the so called halo occ upation 
distribution in some literature (e.g., IPopesso et ai]|2007bl) . is de- 
scribed by a power law, R oc M'\ The correlation of cluster mass 
with the optical luminosity, i.e., the mass-to-light ratio M/L, is 
also described by a power law, M/L oc L" , i.e., M oc L 1+v . In 
FigurefT] we show the correlations between cluster mass and clus- 
ter richness, summed luminosity and GGN/rGGN for 24 nearby 
clusters. The uncertainties of richness R , summed lumin osity L r 
and GGN/rGGN are about 10%-20% dWen et al. 1 120091 . We fit 
the correlations with power-law relations, 

logA/ v i r = (-1.43 ±0.07) + (1.55 ± 0.06) log R, (1) 
logM vir = (-1.77 ±0.08) ± (1.49 ±0.05) log L r , (2) 
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and 

logA'/ vir = (-2.11±0.10) + (2.03±0.08) log(GGiV/r GGJ v).(3) 

Here, M v i r has a unit of 10 14 hT 1 Mq, L r has a unit of 
10 10 h~ 2 Lq. The uncertainty of the estimated cluster mass, 
flogM. is mainly determined by the uncertainties of the intercept 
and the slop e in the logarithm for th ree scaling relations in Equa- 
tion <TTT> — <f3T>- lYee & Ellingsonl J2003h defined B gc to be the ampli- 
tude of galaxy-cluster cross-correlation function and found M v i r oc 
_Bg C 64±0 ' 28 . The slope is in agreement with that of our M v i r to 
GGN/rGGN relation. These scaling relations, Equation <[TJ — <f3j , 
will be used to estimate masses of a complete volume-limited sam- 
ple of clusters in the local universe for cluster mass function. 

We can also use a much larger cluster sample at mediate red- 
shift (~ 0.2 < z < 0.25) for cluster mass function. Some massive 
clusters in this redshift range have their masses estimated (see Ta- 
ble QJ. We obtain masses of 17 clusters in the redshift range of 
0.17 < z < 0.26, of which 10 clusters have more than three esti- 
mates. In Figure[2] we show the correlations between cluster mass 
and cluster richness, summed luminosity and GGN/tggn for the 
17 clusters. Most of them are similarly massive of 10 15 h Mq 
and few have smaller masses, so that it is difficult to determine a 
new scaling relations. Here, we calibrate the mass scaling relations 
by assuming the same slopes of Equation {TJ-([3j and finding the 
offsets. We then get the scaling relations, 

log M vlr = (-1.57 ±0.12) + 1.55 log R, (4) 
log Mvir = (-2.03 ± 0.06) + 1.49 log L r , (5) 
and 

log Mvir = (-2.33 ± 0.11) + 2.031og(GGJV/r GG]v ). (6) 

The uncertainties in Equation {4}-l|6} reflect the scatters of masses 
to the mean relations (dashed line). We notice that the scatter is the 
smallest for the M vlI -L r relation for the high redshift data, because 
clusters with more than three estimates (black dots) are very con- 
sistent with the fitting relation (dashed line). Therefore, the cluster 
masses estimated by the M v i r -L r relation may be more accurate 
than other tracers. The offsets between the relations for samples 
at two redshift ranges may come from the problem of the SDSS 
galaxy data. The sky background level is overestimated for nearby 
bright galaxies (12.5 < r < 15.5), so that galaxies have system- 
aticall y fainter magnitudes by 0.15-0 .2 mag than their true magni- 
tude ( Adelman-McCarthv et al. 2008). This can result in systemat- 
ically lower cluster richness and summed luminosity for clusters of 
0.05 < z < 0.1 than clusters of 0.2 < z < 0.25. The two scaling 
relations are used to samples of clusters at two redshift ranges in- 
dependently. Hence, the systematic bias does not affect the final erg 
values from each sample. 



3 CLUSTER MASS FUNCTIONS 

Assuming a G a ussian distribution of mass fluctuation, 
iPress & Schechterl d 19741) used a linear theory to derive the 
first theoretical expression of cluster mass function, which is in 
agreement with mass functions derived from obser vations and 
nume ri cal simulations within a large mass range (e.g. JWhite et al.l 
1 19931; iReiprich & Bohringed l2002h . Recent simulations show 
slightly more massi ve clusters than the Press & Schechter 
mass function give s ( Shet h & Tormenl 1 19991 ; Ijenkins et al.l 1200 ll : 
IWarren et al1l2006h . In this work, we take the form of the cluster 



mass function as Equation (B4) of ljenkins et~aT1 feOOlh . The mean 
differential comoving number density of dark matter halos is 

dn now P° dlna^ 1 _x n „„-,3.82x ™ 

^ = a316 M^^lrTM- eXp( - [ln,J +0 ' 67] } - (?) 

Here, p = 2.78 x W^Qmh 2 Mq Mpc -3 is the comoving den- 
sity of the universe. M is the halo mass within a radius with a 
mean overdensity of 324 times of the mean density of the universe 
(roughly the virial mass, Mioi, if fi m = 0.3). a 2 (M, z) is the 
variance of the linearly evolved density field smoothed by a spher- 
ical top-hat filter that enclose mass M. Here, <r(M, z) = as x /, 
where as is the present linear rms mass fluctuation on the scale of 
8 h^ 1 Mpc and / is a function of M, z, fl m as well as the Hubble 
constant h, the abundance of baryons fig, and the present cosmic 
microwave background temperature Tomb- dlner -1 /dmM can 
be derived from the expression of a(M, z) (see details of a(M, z) 
in Reiprich & Bohringer 2002). The values of Q m and as are the 
main parameters to define the mass function. The other parameters 
does not strongly affect the results in our analysis, thus can be fixed. 
The as strongly depends on cluster mass function at the high mass 
end. Since the mass function is steep at high mass end, the data 
scatter for mass scaling relations induces more low mass to higher 
mass. Thus, the uncertainty of the mass scaling relation, (Ti ok m, is 
included in the fitting. We re-write the mass function with the un- 
certainty on mass estimate to be the Jenkins function convolved by 
a Gaussian function, 



dn(M) 
dlogM 



dn{M' 
d log M 



\g{logM - logM> logM )dlogM', (8) 



where g(x,a) = e~ x '' '' /2a2 / '(V^tvct). 

First, we use a complete volume-limited sample of rich clus- 
ters (R 16.7, 90% complete) in the local universe (0.05 < z < 
0.1) to determine the cluster mass function. Since the photomet- 
ric redshift was used to identify the cluster member galaxies, the 
absolute magnitudes of member galaxies could have large uncer- 
tainties when the estimated cluster redshift slightly deviates from 
its true redshift. To reduce the uncertainty at low redshift, we use 
the spectroscopic redshifts of clusters if its discriminated members 
are spectroscopically observed. The cluster richness, th e summed 
r-band luminosity and GGN/ tgg n are re-calculated as lWen et al.l 
(2009). In this sample, 56 clusters have richness R 16.7, which 
are used to determine the cluster mass function in the local uni- 
verse. 

We apply the scaling relations of Equation <[TJ — <[3j to these 56 
rich clusters in the local universe and calculate the number of clus- 
ters as a function of mass. Figure[3]shows the cluster mass functions 
and the best fit with Equation ^8). From the probability contours in 
the as-Qm plane for three mass tracers (Figure[4}, we find that the 
as and £l m are coupled in the form of as(£l m /0-3) a = ft. From 
the cluster mass distribution using the mass-richness scaling rela- 
tion, we find 



a» 



03 



0.82 ±0.04. 



(9) 



From the cluster mass distribution using the mass-luminosity scal- 
ing relation, we find 



= 0.90 ±0.04. 



(10) 



From the cluster mass distribution using the mass-GGN/rGGN 
scaling relation, we find 
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Figure 3. Mass function for a sample of 56 rich clusters (R ^ 16.7, 0.05 < z < 0.1). The error bars on the horizontal axis are calculated from the 
uncertainties of Equation (T}-(3), and the error bars on the vertical axis are calculated by Poisson statistics. The solid line is the best fit with the cluster 
mass function of Eq uation (8). The dashed line is the cluster mass function of Equ ation (8) w ith Q m = 0.273 and <rg = 0.813 from the WMAP5 data 
I Koma tsu et al . 2009). Data for the mass functions from Reiprich & Bohringer 12002) and Rin es et al ] <2007t) are plotted for comparison. 
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Figure 4. The probability contour in the ag-Q m plane for three corresponding mass tracers in Figure|3] 68% confidence level for the inner curve and 99% for 
the outer curve. 



OS 



03 



0.83 ±0.04. 



(11) 



During the fitting, we have taken into account only statistical un- 
certainties. Assuming fi m = 0.3, the value of as is 0.82 ± 0.04, 
0.90±0.04 and 0.83±0.04 for masses scaled from cluster richness, 
summed luminosity and GGN/rGGN, respectively. 

We also apply the scaling relations of Equation ©-{6]( to the 
a complete volume-limited sample of 810 rich clusters (R ^ 16.7) 
of 0.2 < z < 0.25 to calculate their masses, and get the cluster 
mass function. Again, spectroscopic redshifts of 466 clusters are 
used since they are available from the SDSS, otherwise photometric 
redshifts are used. Figure [5] shows the cluster mass functions and 



Figure [6] shows the contours in the og-fim plane based on three 
mass tracers. Since there are much more clusters in this sample, the 
mass functions have small errors than those of 0.05 < z < 0.1. We 
fit the data to Equation (8), and find 



03 



03 



= 0.85 ±0.02, 



0.94 ±0.02, 



= 0.82 ±0.02, 



(12) 
(13) 
(14) 
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Figure 5. The same as Figure[3]but for a sample of 810 rich clusters (R ^ 16.7, 0.2 < z < 0.25). The curve from the WMAP5 result by Koma tsuet alj 
1200911 is plotted for comparison. 
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Figure 6. The same as Figure[4]but corresponding to the three mass tracers in Figure[5]for the cluster sample of 0.2 < z < 0.25. 



for the cases using the mass tracer of richness, summed luminosity 
and the GGN/tggn, respectively. Assuming Q m = 0.3, the value 
of a 8 is 0.85 ± 0.02, 0.94 ± 0.02 and 0.82 ± 0.02, respectively. 
They are consistent with those from the cluster sample of 0.05 < 
z < 0.1 for each mass tracer. 



4 DISCUSSIONS AND CONCLUSIONS 

Cluster mass function can be accurately determined from a com- 
plete volume-limited sample. The scaling relations of cluster mass 
have been determined for three optical observations, cluster rich- 
ness, summed luminosity and GGN/rcGN- The scaling relations 
are then used to estimate cluster mass for two samples of rich clus- 
ters. We get cluster mass functions and fit them with a theoretical 



expression. Cosmological parameters are constrained in the form 
of <r 8 (n m /0.3) a = P, with a =0.40-0.50 and f3 =0.8-0.9. For 
fi m = 0.3, we get as =0.8-0.9 using different mass tracers or 
using the rich cluster samples at different redshift ranges. 

The as values from the mass tracers of richness R and 
GGN/tggn obtained using both cluster samples are consistent, 
while as values derived from L r are higher. This discrepancy may 
come from some potential systematic bias on the mass scaling rela- 
tions. If the A/ V i r -L r relations for both samples are really unbiased, 
then the cluster masses tracer by richness R and GGN/rGGN are 
systematically underestimated. However, it is hard to assess which 
one is a better mass tracer. Given the scarce of mass estimates from 
different methods for the same clusters in Table [TJ for the scaling 
relations, it is also hard to estimate the systematic bias on these 
mass estimates due to different methods (X-ray or weak lensing). 
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Table 2. Comparison of results on ag -Q m derived from clust er mass function (upper part) and cosmic microwave background (CMB) mea- 
surement (middle part). See Table 5 of Hetterscheidt et al. (2007) for the results derived from weak leasing statistics. 



Reference Sample No. of clusters crg-f2 m relation <rg (f! m = 0.30) 

or method or observation 



Viana & Liddle (1996) 


X-ray 


25 


n K(i i n 1£;0 n <iko2 
r\ s*nr^~ U.by + U. ±Oii m — U.UDii 

erg = U.oUS2 m 


1.16 


Eke et al. (1996) 


X-ray 


25 


crg=(0.52±0.04)fi~ ' 52+0 ' 13ym 


0.93±0.07 


Markevitch (1998) 


X-ray 


30 


a s = 0.78 ± 0.04 with tt m = 0.30 fixed 


0.78±0.04 


Pen (1998) 


X-ray 


25 


<r 8 =0.53n™ - 53 


1.00 


Borgani et al. (1999) 


X-ray 


70 


<7 8 =0.58 ±0.06n-° 47±016n - 


0.96 


Viana & Liddle (1999) 


X-ray 


10 


<t 8 =0.56Q™°- 47 


0.99 


Blanchard et al. (2000) 


X-ray 


25 


a 8 = 0.96 with O m = 0.30 fixed 


0.96 


Oukbir & Arnaud (2001) 


X-ray 


69 


„ „„-0.57+1.45fi m -3.48n;L+3.77n3 -1.49ni, 

<78 =0.59X! m 


0.91 


VV LI ^ ZUU 1 ) 


X-ray 


~>s 
-j 


a*-0 477rr - 3 -°' 17f2 ™ 34 -°- 13f2 A 


87 


Rnronni pt nl ["90011 
ijUIL^dlll Cl dl. ) 


X-ray 


i n i 

1 UJ 


_ _n fi6+ 06 o - o S5+ - 13 




r lerpaoli et al. ( ZUU 1 ) 


X-ray 


in 

JO 


„ In / iot;+0-034\ n -0.60 
<T 8 =(0.495_ Q37 )S2 m 


i no +0.07 


Viana et al. (2002) 


X-ray 


452 


o-8=0.38S2 m 


0.61 


Reiprich & Bohringer (2002) 


X-ray 


106 


o- g =0.43n™ - 38 


0.68 


Seliak (2002) 


X-ray 


30 


o- g (Q m /0.3)°- 44 = 0.77 ± 0.07 


0.77±0.07 


Viana et al. (2003) 


X-ray 


40 


erg = 0.78+Q 06 with Q ™ = °- 35 fixed 




Schuecker et al. (2003) 


X-ray 


452 


„ -(1711+0 .039 _ n o-i+0.031 
o-g-0.711_ 031 , U m - 0.341_ Q 029 


- 

n 77 +0.05 
u -"-0.04 


Pierpaoli et al. (2003) 


X-ray 


63 


cr 8 =0.77to 04 with O m = 0.30 fixed 


Allen et al. (2003) 


X-ray 


111 


<Tg=(0.508±0.019)n™°' 253±0 ' 024 


0.69±0.05 


Henry (2004) 


X-ray 


51 


erg = 0.66 ± 0.16 with n m = 0.30 fixed 


0.66±0.16 


Dahle (2006) 


X-ray 


35 


CT g(Q m /0.3)°- 37 = 0.67+.™ 4 


n 67+ 04 
u -°'-0.05 


Rines et al. (2007) 


X-ray 


66 


as = 0.92l°; 24 , Q m = 0.24l°;J 4 


0.84 ± 0.03 


Henry et al. (2009) 


X-ray 


48 


cr 8 (!rj m /0.32) - 30 = 0.86 ± 0.04 for C7 m 0.32 
o-g(n m /0.32)°- 41 = 0.86 ± 0.04 for Q m > 0.32 


0.88±0.04 


Vikhlinin et al. (2009) 


X-ray 


49 


<Tg(!TJ m /0.25) ' 47 = 0.813 ± 0.027 


0.75±0.02 


w tiite et al. ( i yy 3 ) 


Optical 




_„ n ^70" 0.56 

as — V.Otilm 


1 1 "> 

1 . 1 L 


Bahcall & Fan (1998) 


Optical 


3 


asn9„ 29 = 0.8 ±0.1 

O 771 


1.13±0.14 


Girardi et al. (1998) 


Optical 


152 


<T8 =(0.60±0.04)fi™°- 46 + 09nm 


1.01±0.07 


Bahcall et al. (2003) 


Optical 


300 


crgfi ^ 60 = 0.33 ± 0.03 


0.68±0.06 


Eke et al. (2006) 


Optical 




n oc:r ,-0.92-4.5(n m -0.22) 2 
crg=0.25Q m 


0.78 


Rozo et al. (2010) 


Optical 


13832 


0"g(f2m/O.25)°' 41 = 0.83 ± 0.03 


0.76±0.03 


Komatsu et al. (2009) 


CMB 


WMAP5 


as = 0.81 ± 0.03, U m = 0.27 ± 0.01 




Liu & Li (2009) 


CMB 


WMAP5 


as = 0.92 ± 0.04, Q m = 0.32 ± 0.03 




Larson et al. (2010) 


CMB 


WMAP7 


a 8 = 0.80 ± 0.03, n m = 0.26 ± 0.01 




Readhead et al. (2004) 


CMB 


CBI 


as = 0.96ig'^ (68% confidence level) 




Dawson et al. (2006) 


CMB 


BIMA 


as = lmtoH (68% confidence level) 




Reichardt et al. (2009) 


CMB 


ACBAR 


„ n qq+0^04 
erg - 0.93I 05 




Sievers et al. (2009) 


CMB 


CBI 


„ n no+0.05 
erg - 0.92I 05 




Savers et al. (2009) 


CMB 


Bolocam 


erg < 1.57 (90% confidence level) 




Veneziani et al. (2009) 


CMB 


BOOMERANG 


erg < 0.92 (95% confidence level) 




This work (0.05 < z < 0.1) 


R 


56 


o-g(T2™/0.3) a42 = 0.82 ± 0.04 


0.82±0.04 








o- 8 (Q m /0.3)°- 46 = 0.90 ± 0.04 


0.90±0.04 




GGN/raoN 




cr 8 (n m /0.3)°- 40 = 0.83 ± 0.04 


0.83±0.04 


This work (0.2 < z < 0.25) 


R 


810 


cr 8 (n m /0.3) - 42 = 0.85 ± 0.02 


0.85±0.02 








o-g(Q m /0.3)°- 46 = 0.94 ± 0.02 


0.94±0.02 




GGN/raGN 




o-g(Q m /0.3) - 39 = 0.82 ± 0.02 


0.82±0.02 



In our work, one potential systematic bias may come from the con- 
version of cluster mass from measured radii to the virial radius. 
Here, we use 7 = M v ir/M v ir,truc to stand for the systematic bias 
of masses in TableQ] where Mvir.true stands for the true virial mass 
of a cluster. Assuming a 7, we get M v j r ,true and then fit the mass 
function of clusters to obtain as- Figure [7] shows the variation of 
erg (with Q. m = 0.3 fixed) as a function 7 based on the M v i T -L r 
relation. We are only concerned about the cases 7^1. For exam- 
ple 7 = 1.3, i.e., masses systematically overestimated by 30%, the 



values of as are lower by about 10%. In fact, the deviation of 7 
from 1.0 is related to the uncertainty of intercept in the logarithm 
scaling relations in Equation ([TJ-((6](. The other possible systematic 
bias on erg may come from the slope uncertainties of the scaling 
relations. Here, we illustrate the dependence of erg on the slope un- 
certainty, Ais. We only apply to the M v i r -L r relation, for example. 
Given a Av, i.e., M vir = A L^+ l/+A ", here v = 1.49 according to 
Equation (2) and ((5}, we fit the power law with the data in Figure[T] 
and [2] to get A, and then get the cluster mass function and fit for 
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Figure 7. The value of erg (with O m = 0.3 fixed) from cluster masses 
based on the Af v ; r -L r relation varies with a possible systematic bias on 
mass conversion of 7 = M v i r /M v i r . true- 
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Figure 8. The value of erg (with Q m = 0.3 fixed) from cluster masses 
based on the M v i T -L r relation varies with a possible systematic bias on 
the slope of the scaling relation by Au. 



erg. Figure [8] shows the erg value varies with Au. We find that the 
erg from the cluster sample of 0.05 < z < 0.1 does not change 
significantly with Au, while the erg decreases from 1.05 to 0.81 for 
the cluster sample of 0.2 < z < 0.25 when the slope varies by Au 
from -0.4 to 0.4. 

We can compare our results of erg with previous determina- 
tions from cluster mass function, as listed in Table [2] Most of pre- 
vious results are based on X-ray flux-limited cluster samples. Our 
results are systematically larger than those from the mass function 
of X -ray cluster s. 

IRozo et al] ( l2010f) used the largest number of clusters from 
SDSS maxBCG catalog to determine the amplitude of cluster mass 
function. They did not estimate the mass for each cluster, but gave 
a statistical mass for clusters within a richness bins by weak lens- 
ing. They got erg = 0.76 ± 0.03 assuming Q. m = 0.30. The 
maxBCG clusters were selected based on the red brightest cluster 
galaxies (BCGs). However, the maxBCG method may miss about 
25 % clusters in which the BCGs have emission line and blue col- 
ors dKoester et al. 2007]). We notice that about 15% rich clusters 
{R ^ 16.7) are missing by the maxBCG method compared to our 
sample in the redshift range of 0.2 < z < 0.25. However, the 



systematic incompleteness only induces an underestimate of 3% 
for erg. Therefore, the discrepancy probably comes from the uncer- 
tainty of mass scaling relations. 

If we take fi,„ = 0.26 derived from WMAP7, then our values 
of erg should become larger by a factor of (0.26/0.3)~ 0,42 = 1.06, 
roughly equal to adding 0.05 to the our erg value in Table|2] There- 
fore, the erg values we derived from galax y clusters are s lightly 
larger than the tho se from the WMAP data jKomatsu et alj|2009l : 
lLarson et a l. 2010). While some reanalysis of the WMAP5 data 
independently iU et al.l |2009|) gives erg = 0.921 ± 0.036 for 
n m = 0.32 ± 0.03 (see lLiu&Lill2009l) . Some studies of cos- 
mic microwave background a t small scales also give higher values 
of erg than that from WMAP dReadhead et alj|2004l ; iDawson et al] 
12004 iReichardt etai]|2009l : Isievers et alj2009h . 

Our result of erg are consistent with many recent studies us- 
ing other methods. For example, the erg by weak lensing method 
ha s a mean value of 0.85± 0.03 (see previous results in Table 5 
of iHetterscheidt et ai1l2007l). which is higher than previous results 
from X-ray clusters. iTegmark et alj d2004l) studied the power spec- 
trum of galaxies from the SDSS to constrain cosmological param- 
eters They obtained erg = 0.89 ± 0.02 and ft m = 0.30 ± 0.03. 
I Led |2009) studied the normalization of the power spectrum via the 
ellipticity function of giant galax y voids fr om SDSS DR5 and ob- 
tained erg = 0.90 ± 0.04. |jena et alj J2005l) used the Lya data and 
found erg = 0.9 and Q m = 0.27. iFeldman et alj j2003h used the 
galaxy peculiar velocities to probe the growth rate of the structure 
and found that erg = 1.13t° j3 ^ fi ™ = O-SOIq ^. 

In this work, we get six values of erg by cluster mass function. 
Basically, the results are consistent. However, the precise value of 
erg is still to be determined since our constraint is not only coupled 
with fl m , but also has large uncertainties on the scaling relations. 
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